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Introduction 


Certes, la croissance de la population Mondiale dans les deux dernières décennies a donné naissance 
à plusieurs problématiques liées à l’urbanisation croissante et à la raréfaction des ressources qui 
induisent le besoin de faire évoluer la gestion des transports ferroviaires urbains en augmentant leur 
efficacité énergétique. Ces problématiques se trouvent à l'intersection des trois piliers du 
développement durable, à savoir l’économie, l'écologie et le social. I| s’agit donc de faire évoluer les 
modes d'exploitation de lignes ferroviaires pour assurer à la fois une amélioration de leur efficacité 
énergétique, une réduction des coûts d’exploitation (ou tout du moins une stabilisation) et une 
maximisation de la satisfaction client. 


C'est donc tout naturellement que j'ai opté le choix pour ce sujet d'optimisation il concerne le 
domaine du Génie Civil. Ensuite il traite une problématique avec des différentes contraintes 
écologiques, économiques, et sociales, au quelles l'ingénieur doit chercher un compromis pour les 
satisfaire. D'autant plus, le Maroc fait partie des pays en cours de croissance, et donc il accorde un 
intérêt particulier au sujet du développement des infrastructures ferroviaires. Ainsi c’est la raison 
principale qui me pousser de se cultiver plus dans l'infrastructure ferroviaire. 


Dans cette Perspective, de Nombreuse recherche ont étés fournîtes afin de proposer des solutions 
pour traiter le problème de l'optimisation énergétique des lignes ferroviaires. Deux grands axes 
d'études sont décrits : la réduction des pics de puissance appelée et la réduction de la 
consommation énergétique globale des lignes ferroviaires. I| en découle de ces deux études que les 
moyens pour réduire la consommation énergétique de lignes de métro sont principalement de trois 
types : la réduction des performances du matériel roulant, l’utilisation de profils d’écoconduite et 
la récupération de l’énergie cinétique issue du freinage 


Dans le présent rapport, je veux me focaliser plus sur la deuxième étude soulignant la nécessité de 
disposer d’un modèle énergétique précis permettant de calculer les flux de puissance qui se 
produisent lors des phases d'exploitation afin de prendre en compte le caractère non linéaire du 
freinage récupératif. Deux résolutions possibles peut-on avoir 


la première : c’est d'optimiser certains paramètres d’exploitations (l'intervalle d'exploitation, les 
temps d’arrêts en station.) sans tenir compte des perturbations du trafic. Ceci en choisissant deux 
méthodes d'optimisation issues de l’Intelligence Artificiel (Optimisation par l’Algorithme Génétique, 
et) 


La deuxième option de résolution est d'apporter une solution au problème de redéfinition en temps 
réel des temps d’arrêt en station pour réduire la consommation énergétique tout en intégrant les 
perturbations de trafic qui se produisent dans les conditions réelles d’exploitation. Ceci en utilisant 
un Algorithme d'apprentissage par renforcement avec un réseau de neurones. 


Par la suite j'en parlerais d’autres solutions qui peuvent être sujet d’élargissement l'efficacité 
énergétique des lignes ferroviaires à savoir l'installation de systèmes de stockage fixes ou 
embarqués ou l'installation de sous-stations réversibles. 


Et avant de terminer, j'indique les opportunités et les enjeux peuvent avoir lieux lors de l'application 
de ce projet d'optimisation au Maroc 


Finalement je termine par une conclusion sur le projet 


Modélistion du Problème : 


1. Fonctionnement d’une ligne de Métro : 


L'optimisation énergétique d’une ligne de métro automatique, nécessite d’abord de 
comprendre comment les rames de métro sont alimentées en énergie et comment les 
échanges de puissances s'effectuent entre les trains et les sous-stations d'alimentation. 

Les principaux éléments qui constituent une ligne de métro automatique sont passés en revue : le 
matériel roulant, les barres de guidage qui permettent d'alimenter les rames en énergie et les sous- 
stations d’alimentation 





1.1 Présentation du réseau de traction (sous station d’alimentation) : 

La figure ci-dessous, présente le schéma d'alimentation électrique d’une ligne de métro, ainsi que les 
niveaux de tension vus par les différents équipements. Les Postes de Livraisons (Plv) distribuent 
l'énergie électrique aux Postes de Redressement (PR) qui alimentent la ligne de métro et aux Postes 
Éclairage et Force (PEF) qui alimentent les stations 


Réseau de distribution HTA 


20kV AC 


20kV AC Poste de Livraison (PLv) 20kV AC 
Poste de Redressement (PKR Post l t Fo PEF 
ay pe [Poste Eclairage et Force (PEF)| 


| 750V DC 


Barres de guidage 


220V AC | 








Schéma de l'alimentation générale d’une ligne de métro 


1.2 Présentation du matériel roulant : 


Schéma d'alimentation des moteurs d’un véhicule. La figure ci-dessous représente le schéma 
électrique simplifié du système de propulsion et freinage d’un véhicule (PBS pour Propulsion and 
Braking System), sur lequel sont indiqués les principaux éléments permettant l’alimentation en 
courant des moteurs synchrones. 







Coffret 
Haute 
Tension 


Schéma d'alimentation des moteurs d’un véhicule 


2. Modélisation énergétique du matériel roulant : 


2.1 Modélisation du déplacement des trains 

Selon le niveau de précision souhaité, il existe classiquement deux approches pour modéliser 
le déplacement des trains sur une ligne ferroviaire : l'approche temporelle et l’approche 
épisodique (le terme épisodique provient de la traduction du terme anglais event-based) 


L'approche temporelle : définit le mouvement des trains par une séquence d'événements 
caractéristiques tels que l’arrivée et le départ de station. L'intérêt majeur de cette approche 
est de pouvoir avancer l'horloge de simulation jusqu’à l'occurrence de l'événement suivant. 


L'approche épisodique : L'approche temporelle consiste à discrétiser l'horizon de 
simulation afin d'étudier le mouvement des trains à chaque pas de temps. 


V(t + St) = V (E) + 7(6).6t 
(X(t + ôt) = X(t) + V(t).ót 


2.2 Modélisation mécanique du matériel roulant 

La modélisation des efforts de traction se fait quasi-systématiquement par application de la 
deuxième loi du mouvement de Newton, appelée aussi Principe Fondamental de la 
Dynamique (PFD) 

| est nécessaire d'identifier les différentes forces qui s'opposent à l'avancement d’un train. 
Ces efforts sont principalement de trois types : les efforts liés aux frottements (contact 


roues-rails et frottements aérodynamiques), les efforts liés à la pente de la voie et les 
efforts liés à l’accélération du train. 
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Profil d'accélération pour un parcours inter station (vitesse 
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Les forces de frottements sont déterminées par l'équation empirique : 


Frsl) = A+ B - V(t)+C - (VE) - W(#)}(#) 


À : représentent respectivement le coefficient de frottements secs (proportionnel à la charge 
par essieu, en N) 

B : le coefficient de frottements visqueux (en N.s.m'!) 

C : le coefficient de résistance aérodynamique (en N.s? .m*°) 

V : la vitesse du véhicule 

W : la vitesse du vent (en m.st) 


Les efforts dus à la déclivité sont exprimés par la relation 


Frente(t) = M.g.a(t) 


Où a, M et g correspondent respectivement au profil topographique de la ligne, à la masse 
du train passagers inclus et à l'accélération de la pesanteur 
Les efforts dus à l'accélération sont quant à eux exprimés par la relation 


Pi (£) = Minertie- M. (t) 


Le coefficient n permet de prendre en compte l'inertie des masses tournantes tandis que y 
désigne l'accélération du train considéré. 


2.3 Modélisation électrique du matériel roulant 


Fetee(t) = n" > V(t) > (Face(t) + Fpente(t) + Frea(t)) + Faux 
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3. Modélisation d’une ligne de métro : 


Le schéma électrique d’une ligne de métro simplifiée est présenté à la figure 2.17. Cette ligne fictive 
est composée de 3 sous-stations d'alimentation (nœuds 1,3 et 5), d’un train en traction (nœud 2) et 
d’un autre en freinage (nœud 4). Le train freineur est distinguable du fait de la résistance variable 


placée en parallèle de la source de courant afin de modéliser la possibilité de dissipation du freinage 
électrique 


3.1 Application à un exemple simplifié 
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3.2 Analyse nodale modifiée 


L'analyse nodale modifiée repose sur l’utilisation des lois de Kirchhoff pour déterminer les 
tensions nodales et les courants circulant dans les branches d’un circuit électrique 


Et En combinant les lois de Kirchhoff et la loi d'Ohm, qui relient l'intensité du courant 
électrique traversant un dipôle avec la tension à ses bornes, il est possible d'écrire les 


équations reliant les éléments électriques entre eux. 


Nous obtenons alors la formulation matricielle suivante où [V] et [I] sont respectivement les 
vecteurs des tensions et des courants nodaux et [Y] est la matrice d’admittance qui synthétise 


la connaissance de l’état du réseau 


FE 


où 


[Y).V] = (A 
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4. Résolution d’un problème de répartition des charges : 


4.1 Formulation mathématique du problème 


L'objectif de la résolution est de déterminer les solutions de l'équation du réseau où Vn et In sont 
respectivement les vecteurs de tensions et de courants aux nœuds du réseau et Yn est la matrice 
d’admittance du réseau. La matrice jacobienne se déduit alors immédiatement 


F(V) =Y, V, -I =0 


F'(V) = Yn 
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4.2 Algorithme de résolution : Méthode Newton-Raphson 
L’algorithme de résolution est synthétisé dans le synoptique de l'algorithme 1. Il s’agit d’un 
processus itératif composé de 3 étapes : le calcul de la matrice jacobienne de la fonction à 
annuler, l'inversion de cette matrice puis le calcul d’une nouvelle solution. Cette relative 
simplicité d’implémentation est également l’un des facteurs ayant contribué à la 
popularisation de la méthode de Newton-Raphson. 


Algorithme 1 Algorithme de Newton-KRaphson 


Création d'une solution initiale V{!) 
Tant que |V"+! — V”| > e Faire 
Calcul de la matrice jacobienne Jp 
Inversion de la matrice jacobienne Jp 
Calcul de la nouvelle solution V” = V-t! — J-F (Va-1) 


Fin du Tant que 


4.3 Résolution par heuristique itérative : 
La méthodologie présentée dans ces travaux a donc été reprise pour pouvoir com parer les 
performances de cette heuristique avec la méthode de Newton-Raph 


à r 
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En outre, le déplacement de chaque train impose des conditions de charge, ce qui a pour effet de 
conditionner les échanges de puissance qui s'effectuent entre les sous stations et les trains 


, 7 
Firas — Virainas-ltrains 


La détermination des paramètres électriques en chaque nœud du réseau se fait ensuite par un 
procédé itératif, résumé par l'algorithme 2 


Algorithme 2 Algorithme d'heuristique itérative 


1: Création d'une solution initiale V® 


2 Tant que [V"*!- V”| > e Faire 


3 Calcul du vecteur courant 7, = Y.V, 
d: | Eaa LE Ín 

| a F | 
5: Calcul du nouveau vecteur tension Virainas = LL. L 
` | a | p Paimas 
ü: |} fl — "rains 
7: Fin du Tant que 


Optimisation des Paramètres d'exploitation 


1. Formulation du problème d’optimisation : 


Objectifs è Optimisation énergétique de tables horaires pré-existantes 
èe Nombre de trains en ligne 





ne è Marge de variation des temps d'arrêt en station 

Moyens è Sélection de l'intervalle d'exploitation le plus favorable 

d'action è Sélection de la combinaison de temps d'arrêt en station la plus 
favorable 

Indicateurs è Taux de réutilisation du freinage électrique 


1.2 Définition des variables utilisées : 
Chaque train du carrousel parcourt les N station stations de la ligne de métro, en suivant des 
profils de vitesse entre chaque station et en effectuant des arrêts en station d’une durée si,) 
(où i = [1, ..., N station] et j = [1, ..., N trains]). Pour les besoins de l'étude, nous considérons 
un domaine temporel discret avec un pas d’échantillonnage variable selon le niveau de détail 
requis. 


1.3 Définition des contraintes : 


Le temps de parcours TP est calculé comme la somme des temps de parcours inter station ti 
additionnée de la somme des temps d'arrêt en station sij et du temps de battement tb) 
IN atatior 


Tp - = > FiF 4 Èi LE E EL, Niraine |]; WE = I 2. N station | 


t— 1 


Le temps de parcours est contraint de sorte que les trains respectent la table horaire nominale pour 
assurer la qualité de service souhaitée par l'exploitant. La durée des parcours inter station étant fixe, 
il est donc nécessaire d'imposer une contrainte sur la plage de variation des temps d’arrêt en station 
ou si, nôm représente le temps d'arrêt nominal que doit effectuer un train à la station i, tandis que 
les paramètres As1i et As2i permettent de modifier l'amplitude de la modulation 


Si nom AS]; < Sij < Si nom ` ÅS; 


l'intervalle d'exploitation I est également contraint pour assurer les exigences en matière de sécurité 
et de qualité de service 





Lo <I< Ip 
INtrains Nirains — 1 


Le paramètre As2i devra alors respecter la contrainte afin que la contrainte (3.3) soit toujours 
respectée 


Is alto 
D (Asa) < Tp — T * Nirains 


i=l 


1.4 Définition de la fonction objective : 


L'énergie consommée par le carrousel sur la durée de simulation est alors définie par : 


Ntroins TP . 
Palt) = $} P(t) ia à 7 


j=l 
seul le critère portant sur la consommation globale d'énergie Etot est étudié. La fonction objective 
est alors donnée par (3.9) et a pour unique but de minimiser la consommation globale d'énergie. 


Fobjectis = min Ftalt) 
OLÉSTy 


2. Optimisation de l'intervalle : 


La première étape dans la conception de tables horaires consiste donc à déterminer pour 
chaque carrousel de trains, l'intervalle d'exploitation permettant de minimiser la 


consommation énergétique 
Algorithme 3 Simulation du parcours type d'un train sur un tour de boucle 


Entrée(s) temps.arret et profilinterstation 
i: Pour i = 1 à Natation Faire 
2 profil = | Piur © longueur(temps.arret(i)): profilinterstation(i) | 
3 parcours.type = | parcours. type; profil ] 
4 Fin du Pour 
Sortie(s) parcours.type 
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Position (m) 


18 2 22 24 26 


1,4 1.6 





0.8 1 1.2 
Durée d'exploitation (8) 


Évolution de position d’un train sur un tour de boucle 


0 02 04 0,6 


Energie dissipée (kWh) 


3. Étude des intervalles d'exploitation : 
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4. Optimisation des temps d’arrêt en station : 


L'objectif de la modulation des temps d’arrêt en station est de modifier la longueur des 
phases d’arrêt en station (phases rouges) pour faire coïncider le maximum de phases 
d'accélération (phases bleues) avec des phases de freinage (phases vertes) afin de maximiser 
la réutilisation de l'énergie électrique générée lors du freinage. Les profils présentés en 
figure 3.10 sont des profils théoriques puisque les phases de freinage sont constituées de 
paraboles indiquant un renvoi total de l'énergie issue du freinage ; dans un cas réel 
d'exploitation, l'allure de la puissance électrique renvoyée est généralement beaucoup plus 
irrégulière puisque les conditions de réceptivité totale de la ligne sont rarement rencontrées. 
Dans le cadre de l'optimisation des temps d'arrêt en station, une solution est définie comme 
la combinaison des temps d’arrêt effectués par les trains en exploitation. Le nombre de 





Nombre de trains en ligne 


solution possibles est déterminé par (3.12) où Amoa représente l'amplitude de la modulation 
et Nstop le nombre d'arrêt en station à optimiser 
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4.1 Détermination de la méthode dď’optimisation : 


LI 
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4.1.1 Optimisation par métaheuristique 


L’optimisation s'effectue par une succession d’itérations appelées générations et consiste à faire 
évoluer une population d'individus grâce à des opérateurs de sélection naturelle pour générer une 


nouvelle génération d'individus. Les performances des individus sont alors évaluées par la valeur de 
leur fonction objectif ou fitness qui reflète la capacité des individus à voir leurs génotypes être 
propagés à la génération suivante. Plus généralement, la fonction objective caractérise le degré 
d'atteinte de l'objectif de l'optimisation. 


4.1.2 Optimisation par algorithme génétique 


Population initiale de Npop individus 


Evaluation des individus 


Opérateur de sélection 





Opérateur de croisement 





Opérateur de mutation 





Insertion de nouveaux individus 


Algorithme 5 Algorithme d'optimisation par algorithme génétique 
Entrée(s) Initialiser une population de Npop individus 
1: Tant que Le critère d'arrêt n'est pas vérifié Faire 
2 Evaluer la valeur de fonction objectif de chacun des Npop Individus de la popu- 
lation. 
3: Appliquer un opérateur de sélection pour choisir les Npar parents de la génération 
suivante 
4 Sélectionner des couples de parents et appliquer un opérateur de croisement avec 
une probabilité Fas, pour générer des couples d'enfants 
&  Apphquer un opérateur de mutation à chacun des enfants avec une probabilité 
Fon 
6 Fin du Tant que 
Sortie(s) Population finale Npop individus 


Algorithme 6 Sélection "Roulette" 
Npop | 
l: S1 = > Fevai(Ci) p 
2 Choisir un nombre L tel que L € [0; 51] . 
3: Tant que $: < L Faire 
4 1 1+1 
Sat Sa + Fival(Ci) 
. Fin du Tant que 


. Return t 


=] oh 


4.1.3 Optimisation par algorithme génétique 


Dans la pratique, les particules de l’essaim sont considérées comme des solutions possibles au 
problème. Une particule est définie par une position Xn i sur l’espace des solutions et une vitesse de 
déplacement V n i 3 , elle possède également une mémoire qui stocke la meilleure position visitée 
Pi,best et la valeur de la fonction objectif en ce point. Les particules ont également la capacité de se 
communiquer entre elles la position de la meilleure solution globale connue de l’essaim Gbest 


L’équation du mouvement est donnée par le système (4.30). wv n i représente l’inertie de la 
particule, r1B1(Pi,best - Xn i ) est la mémoire cognitive de la particule, tandis que le dernier terme 
r2B2(Gbest - Xn i ) exprime la mémoire sociale de l’essaim. Le coefficient d'inertie w sert à contrôler 
l'importance de la direction de recherche courante sur le déplacement futur et pour régler la 
capacité d'exploration de l’essaim tandis que B1 et 62, respectivement le coefficient cognitif et le 
coefficient social, permettent de régler la capacité d'exploitation. r1 et r2 sont quant à eux des réels 
tirés uniformément dans l'intervalle [0, 1]. L'ensemble de ces coefficients a donc pour objectif 
d'orienter la recherche de l’optimum 
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Algorithme 7 Algorithme d'optimisation par essaims particulaires 
1: Imitialiser une population de Wpop individus, définis par des vitesses aléatoires et 
des positions prises dans l'espace des solutions. 
2: Evaluer les valeurs de fonction objectif aux positions des Npop individus. 
3: Définir les optima initiaux Fi beat €t Cheat. 
4: Tant que Le critère d'arrêt n'est pas vérifié Faire 


5: Pour i= 1 à Npop Faire 

6: Mettre à jour la vitesse de la particule suivant l'équation (3.144) 

7: Mettre à jour la position de la particule suivant l'équation (3.14b) 

8: Evaluer la valeur de la fonction objectif à la position de la particule f(X;(n)) 
9: Si f(Xi(n)) < J (Fimea) Alors 

10: Pi beat + Xlr) 

11: Fin du Si 

12: Si f(Xi(n)) < fGen) Alors 

13: Gbe + Xilin) 

14: Fin du Si 


15: Fin du Pour 
16: Fin du Tant que 


4.1.4 Hybridation des méthodes d'optimisation : 


Le principe d’hybridation retenu est qu’à chaque itération, la population issue de l’optimisation par 
OEP a connaissance de la meilleure position visitée par la population de l'AG, et inversement, l'AG a 
la capacité d'utiliser les solutions déterminées par OEP pour former de nouveaux individus. De cette 
manière, chaque population bénéficie des propriétés de recherche et de convergence de l’autre 
méthode. 


Cette hybridation a été privilégiée pour plusieurs raisons. D’une part, l’utilisation de deux 
populations générées par deux méthodes d'optimisation différentes permet de conserver les 
propriétés de convergence de chacune des méthodes à chaque itération, d'autre part, le partage des 
informations entre OEP et AG augmente la vitesse de convergence globale et la diversité des 
solutions explorées. 


Initialisation des paramètres 





Création de la population initiale 





Evaluation des individus par MON 










Déterminer F; peut Et Gres 


mon — -+ = - 
Classer les individus par fitness 


Opérateur de sélection 
Mettre à jour 
les vitesses = | = 
. Opérateur de croisement 
Eet FOsitlons 
des individus 


Opérateur de mutation 


Insertion de nouveaux individus 


A — _ — —— _____ Pac 
Fusion des populations Forer et Par 


[Foer Fac] 


5. Comparaison des méthodes d’optimisation : 


Chacune des métaheuristiques présentées précédemment : OEP, AG et l’hybride OEP-AG a 
été implémentée en vue de solutionner le problème d'optimisation des temps d’arrêt en 
station. Ces méthodes ont ensuite été comparées pour déterminer laquelle est la plus 
efficace pour le type de problème rencontré 


L'indicateur de performance utilisé pour comparer ces techniques est le gain énergétique 
entre une table horaire utilisant des temps d’arrêt en station nominaux et le meilleur 
planning de stationnement trouvé par optimisation. 

Le gain énergétique est donc défini par rapport à la consommation du carrousel sur un tour 
de boucle en considérant un fonctionnement nominal 

Le tableau ci-dessous présente une comparaison des performances de convergence des 
populations de [OEP et de l'AG et la densité de distribution des solutions renvoyées par 
chaque méthode. La colonne de gauche reprend le cas où chaque population évolue 
indépendamment l’une de l’autre ; tandis que la colonne du milieu correspond à l’évolution 
des populations dans le cadre de l’hybridation 
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Gain énergétique (9) 


| Méthode AG | Méthode OEP | Méthode hybride OEP-AG 
Temps de calcul - 


unitaire (s) 


G.8613 6.7158 ERESI 





Élargissement Du sujet 


1. Optimisation temps réel des tables journalière en 
tenant compte des aléas du Traffic 


Les techniques d'optimisation par AG et OEP ont alors été développées pour assurer la 
résolution du problème, cependant, la nécessité d'assurer le compromis exploration 
exploitation en un temps de calcul le plus faible possible a amené à définir une méthode 
hybride permettant de synthétiser les points forts de AG et de l’OEP dans un seul 
algorithme de recherche. Moyennant une détérioration des performances en temps de 
calcul, cette méthode hybride permet alors d'augmenter grandement la densité de 
solutions jugées acceptables tout en assurant une convergence vers un optimum plus 
élevé que ceux trouvés par les méthodes AG et OEP 

Cependant, il est encore nécessaire d'apporter une solution au problème de 

re définition en temps réel des temps d’arrêt en station pour réduire la consommation 
énergétique tout en intégrant les perturbations de trafic qui se produisent dans les 
conditions réelles d'exploitation. 


2. Stockage énergétique fixe ou embarqué 


Ils existent pour améliorer encore l'efficacité énergétique des lignes ferroviaires comme 
l'installation de systèmes de stockage fixes ou embarqués ou l'installation de sous- 
stations réversibles. 

Le stockage énergétique stationnaire consiste en une ou plusieurs unités de systèmes de 
stockage placés le long de la ligne, dont le but est de récupérer l'énergie électrique issue 
du freinage générée en excès sur la ligne [32]. Le stockage énergétique embarqué 
consiste quant à lui à équiper les trains en exploitation de systèmes de stockage. Cela 
permet aux trains de stocker l'énergie issue du freinage qui n’a pas pu être consommée 
par les trains présents sur la ligne. Cette dernière solution présente l'intérêt majeur de 
ne pas impliquer de pertes en lignes supplémentaires liées au renvoi de l'énergie du 
freinage entre les trains, mais également d'assurer une certaine autonomie énergétique 
des trains. 


Applicabilité Au Maroc 


Certes le Maroc, a accorder une importance particulière au sujet du développement des 
infrastructures ferroviaire au cours de la dernière décennie. 


D'ailleurs ces travaux d'optimisation sont focalisés sur la conception d’une méthode logicielle 
permettant d'améliorer l'efficacité énergétique de lignes ferroviaires. Cette solution présente 
l'avantage de ne pas nécessiter d'investissement lourd. Ainsi il s’agit d’une grande opportunité 
d'appliquer ces techniques d'optimisation au Maroc pour enrichissez dans un premier temps le 
système de transport dans les grands pôles en préservant gain énergétique énorme 


Et dans un deuxième temps de contribuer avec plusieurs organismes spécialisés dans le domaine à 
savoir les bureaux d'étude pour que ces projets avaient plus d'efficacité et de fiabilité scientifique 


Conclusion 


A travers cette recherche, Nous avons détaillé Les différentes méthodes d'optimisation permettant 
d'implémenter un outil d'aide à la décision permettant de resynchroniser en temps réel les trains de 
manière à limiter la dissipation de l'énergie issue du freinage. 

Un modèle énergétique d’une ligne de métro a d’abord été défini, puis une méthode de résolution 
itérative a été décrite pour pouvoir calculer les flux de puissances effectués entre les trains et les sous- 
stations. Ensuite, une méthode d'optimisation hybride a été développée pour déterminer des 
combinaisons de temps d’arrêt en station permettant de minimiser la consommation énergétique. 
Puis, la méthode de résolution itérative a été synthétisée par un réseau de neurones afin de gagner en 
temps de calcul tout en conservant la précision sur les flux de puissances réels. 

D'ailleurs, Les développements futurs de ces travaux de thèse pourraient consister d’une part à 
développer une modélisation de lignes de métro encore plus précise afin de se passer de la nécessité 
de posséder des enregistrements réels pour en déduire une politique optimale pertinente, et d'autre 
part à développer une méthode d'optimisation mathématique des temps de stationnement pour 
s'assurer de constituer une base de données présentant des solutions optimales. 
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